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Motivated by recent studies of inclusive neutrino nucleus processes and muon capture within a 
correlated local Fermi gas model (LFG), we discuss the relevance of nuclear finite size effects in 
these reactions at low energy, in particular for muon capture. To disentangle these effects from 
others coming from the reaction dynamics we employ here a simple uncorrelated shell model that 
embodies the typical finite size content of the problem. The integrated decay widths of muon 
atoms calculated with this shell model are then compared for several nuclei with those obtained 
within the uncorrelated LFG, using in both models exactly the same theoretical ingredients and 
parameters. We find that the two predictions are in quite good agreement, within 1-7%, when the 
shell model density and the correct energy balance is used as input in the LFG calculation. The 
present study indicates that, despite the low excitation energies involved in the reaction, integrated 
inclusive observables, like the total muon capture width, are quite independent of the fine details of 
the nuclear wave functions. 

PACS numbers: 23.40.Bw: 25.30.-c; 21.60.Cs; 



I. INTRODUCTION 

In this paper we study the importance of nuclear fi- 
nite size effects in inclusive muon capture reactions. The 
motivation for this investigation comes from the results 
of a recent publication 0, 0, where we have devel- 
oped a model which describes rather well the inclusive 
12 C(z^,/i~) and C(y e ,e~) cross sections near thresh- 
old, and inclusive muon capture by nuclei. This ap- 
proach, which is an extension of the quasi-elastic inclusive 
electron scattering model of , is based on a Local Fermi 
Gas (LFG), where the simplicity of the model makes it 
possible to include agreat variety of effects into the re- 
action dynamics 0, UIE 0, ■ In particular long-range 
nuclear correlations are taken into account by comput- 
ing the RPA polarization propagator containing nuclc- 
onic and A(1232) degrees of freedom. All these effects 
are crucial for the correct analysis of atmospheric neu- 
trino fluxes la. flCJ. and to describe the recent neutrino 
experiments gUSlll EE EE El 

The results of pj , particularly those of muon capture, 
indicate that for some kind of inclusive reactions the va- 
lidity of the LFG can be extended to energies lower than 
expected. Although the LFG leads to reasonable predic- 
tions for integrated quantities, at low energies it is not 
possible to describe the shape of the differential neutrino 
cross section or muon capture width within this model. 
In fact both the discrete and the continuum states of 
the final nucleus, including giant resonances, contribute. 
However, when one sums over all the final states, the in- 
formation about the fine details of the spectrum is lost 
and only the global contribution remains. Therefore for 
some integrated inclusive observables, the results depend 
mainly on global quantities such as the correct energy 
balance or the nucleon distribution. An example is the 
inclusive pion capture model of 0, Il9|. I t is also worth 
to mention the pioneering work of ref. |20| in the context 



of comparing shell model with Fermi gas for u-inclusive 
neutrino scattering. 

The goal of this paper is to investigate whether finite 
nucleus effects can affect significantly the LFG results of 
ref. for inclusive muon capture. From the present 
calculation it is possible to estimate an uncertainty of 1- 
7% due to finite size effects not taken into account in the 
LFG calculation of There already exist microscopic 
calculations of neutrino-nucleus reactions and muon cap- 
ture, based on the RPA or large shell model (SM) basis 
ISmiHIHIil El, and other approaches such as 
the relativistic shell model j2^ and the Green's function 
method |2^|. All of them treat correctly the finite size 
of the system. However from direct comparison of these 
models with the LFG it is not possible to deduce the 
importance of the finite nucleus treatment. The reason 
is that such comparison should be done between mod- 
els embodying exactly the same nuclear dynamics. In- 
stead a great variety of residual interactions, shell model 
wave functions, current operators, Coulomb effects, etc., 
have been used in these works, making it impossible to 
disentangle the impact of the finite size on the different 
contributions. 

Before such comparison be undertaken in a sophis- 
ticated framework, it is convenient to understand the 
uncorrelated case. For this reason we have chosen to 
perform this comparison with a simple model where the 
finite-size effects can be easily recognized. Thus here we 
consider the extreme SM, i.e., single particle states in a 
Woods-Saxon (WS) potential, and we compare the re- 
sults with those obtained from a LFG model. We do 
not include long range correlations of RPA type or con- 
figuration mixing. In addition we use the static form 
of the single- nucleon charged current (CC), in order to 
simplify even more the calculation and to draw cleaner 
conclusions. On the other hand, within our simplified 
approach, we do not confine ourselves to the single case 
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of C, but make a more thorough study of how finite size 
effects vary with increasing nuclear mass, by comparing 
the results for a set of closed shell nuclei: 12 C, 16 0, 40 Ca, 
and 208 Pb. 

The present calculation represents a first-stage test 
case to check the "equivalence" of LFG and shell models 
for some inclusive processes. We choose the /i-capture 
reaction for this investigation since it involves low ex- 
citation energies, the worst conditions for the LFG. Of 
course, under these simplifications it makes no sense to 
compare our results with the experimental data |30| , nor 
it is the intention of this work, since it was already done 
in ref. with the full model. 



II. GENERAL FORMALISM 
A. Partial and differential width 

Here we present the formalism to describe the inclusive 
muon capture within our model. We use Bjorken& Drell 
[3l| conventions. We consider a negative muon bound 
into an initial nucleus ^X, which decay into a final nu- 
cleus z~iY pl us a m uon neutrino (not detected). 
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The final nucleus can be in the discrete or in the con- 
tinuum. We assume that the initial muon is in a s-wave 
state (normalized to one) 



(jiis (r) = 4> ls (r) 



Risir) 



47T 



(2) 



We describe the wave function <pi s in a non-relativistic 
framework by solving the Schrodinger equation for the 
muon in the nuclear Coulomb potential, including finite 
size and vacuum polarization effects. The final neutrino 
has four- momentum fc' M = (e',k'). The leptonic current 
matrix element involved in the decay is then 



(^\f(x)\fi) =£^ ls (x)e 



ik x 



(3) 



where x^ — (t,r) is the space-time coordinate, <j)i s (x) is 
the time-depending muon wave function 



4>u(x) = 4>\ s (r)e~ 



(4) 



and e is its initial energy (including the binding). Since 
we treat the muon as non-relativistic, we describe its spin 
by a Pauli spinor \ which is contained into the leptonic 
vector defined by 



Ve' 



1/2 



^(k')7 A '(l-75)^(0) 



(5) 



where we have written the muon four-spinor as u^O) = 
(x, 0), i.e., corresponding effectively to a four-spinor with 



momentum zero. This is equivalent to neglect in the fol- 
lowing the initial muon momentum k = in the kine- 
matics — however the full spatial dependence of the wave 
function 4>i s {r) is maintained in the matrix element, see 
below. Finally in JSJ) V is the normalization volume of 
the neutrino plane wave, and m! its mass, that is set to 
zero at the end of the calculation. 

The S'-matrix element relevant for the decay reaction 
JTJ) is then 



S fi = -2iri5(E f - E, 



W )^(/|J„(- 



-kOli) (6) 



where \i) and |/) are the initial and final nuclear states, 
with energies -E 1 , and Ef respectively, u> = e— e! is the en- 
ergy transfer, G — 1.1 664 x f _5 GeV -2 cos C is the Fermi 
coupling constant multiplied by the cosine of Cabibbo's 
angle, and we have introduced the effective current oper- 
ator J M , defined in coordinate space as 



^(r) = J M (r)0i s (r). 



(7) 



Here J M (r) is the nuclear CC operator to be specified 
below. Finally, J M (q) is the Fourier transform 



J M (q) = / d 3 re^ r J M (r). 



(8) 



The differential decay width can be computed easily. 
Since the initial nucleus is unpolarized, the distribution 
of neutrinos is independent of the angles, and the corre- 
sponding angular integral gives a factor 4ir. 

At this point we have to distinguish two cases, depend- 
ing on the kind of final state |/) reached. In the first case 
the final nucleus is in a discrete state, that can be the 
ground state or an excited state. The neutrino energy 
takes discrete values fixed by energy conservation, and 
the partial width for the transition i — > / is written as, 



G 2 e' 
27r m 



(9) 



where q = 
introduced 



jq|, and the usual leptonic tensor has been 



jfw = k n k >» + k v k <n _ me ' g ^ + i e ^ kak ' (io) 

for initial muon momentum — (m, 0), where m is the 
muon mass. We have also defined the muon-hadronic 
tensor for the transition 



wfc f (q)= J2 </l^(q)K>*</l-Mq)l*> 

M f Mi 



(11) 



where q = k — k' = — k' is the momentum transfer, we 
sum over final spin components My, and average over 
initial spins Mi. 

In the second case, the final nucleus goes to the contin- 
uum, above the one-particle emission threshold, and the 
final neutrino energy ranges between and the maximum 
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energy available minus the nucleon separation energy of 
the final nucleus. The continuum spectrum of neutrinos 
is described by the differential decay width 



dT c G 2 e' 



r,» v W$(q,u). 



(12) 



de' 2tt m 

where now the continuum hadronic tensor is defined as 
W${q,w) = Y^5{E f - Ei w)(/|J p (q)|*)*(/|J,(q)|<). 

fi 

(13) 

Here a sum over final (continuum) states and an average 
over initial spin is assumed. 

The contraction between the leptonic and muon- 
hadronic tensor is easily performed in a coordinate sys- 
tem where the z axis is in the q direction. We finally 
obtain the following expression for the differential decay 
width 



dT c _ G 2 
~dV ~ 2^ 



e' 2 (R c + R L - 2Rcl + R T + 2Rt> ) (14) 



and a similar expression for the discrete partial widths, 
where for simplicity the response functions have been in- 
troduced as the following components of the hadronic 
tensor HI HI: 



Rc 


= w 00 


(15) 


Rcl 


= -1(^3 + ^3°) 


(16) 


Rl 


= w 33 


(17) 


Rt 


= w n + w 22 


(18) 


Rt> 


= - l -(w 12 -w 21 ) 


(19) 



The total (inclusive) width is obtained by integrating 
and summing over the continuum and discrete, respec- 
tively 



r = £r w+ f 



dTc 
de' 



de' 



(20) 



B. Multipole expansion 

Since the shell model states have good angular momen- 
tum, \i) = \JiMi), \f) = \JfMf), it is usual to perform 
analytically the sums over third components using the 
Wigner-Eckart theorem. To this end one begins with the 
following multipole expansion valid for the components 
of any current operator in momentum space as a sum of 
operators with good angular momentum of rank J (note 
that the z-axis is in the q direction) 



J (q) = V^J2i J [J]C J0 (q) 



(21) 



J m {q) 



4^]>> 7 [J]L J0 ( 9 ) 

oo 



J=0 



Ej m + mM Jm (q) 



(22) 



±1 



(23) 

where we use the notation [J] = y2J + 1, and in the last 
equation the spherical components of the current vec- 
tor have been introduced J±i = t(Jx ± J y )/*j2. The 
operators in this expansion are the usual Coulomb, lon- 
gitudinal, transverse electric and transverse magnetic op- 
erators, defined by 



Cj (q) = / d 3 rjj(qr)Yjo(r)J (r) 



(24) 



Ljo(q) = ~ / d 3 rV[jj(qr)Yj (r)}-J(r) (25) 

q J 

Ej m (q) = - f d 3 rV x \jj(qr)Yj Jm (r)}-3(r)(26) 

q J 



M Jm (q) = d 3 rjj(qr)Yj Jm (v)-~3(r) 



(27) 



where jj is a spherical Bessel Function and Yjj m is a 
vector spherical harmonic. Note that the above expan- 
sions H21H23fl are a direct consequence of the familiar 
plane wave expansion in spherical Bessel functions and 
spherical harmonics, inside the Fourier transform JSJ. 

Inserting the expansions H21H23J1 inside the hadronic 
tensor (|13f) we obtain 

Rc = ^^T.\Cj\ 2 (28) 



2J t 



Rl 



Rcl 



Rt 



4tt 



2J 4 4 
2n 



^>L.j\ 2 



2J % 



1 ^ 

J 



(C*jLj + L*jCj) 
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23, 



1 ^ 

J 



\Ej\ 2 + \Mj\ 2 ) 



Rt> — — 



2tt 



2J, 



(29) 
(30) 

(31) 
(32) 



for the responses in the discrete, and a similar expression 
for the continuum responses with the addition of a sum 
over final states and a delta of energies J^f 5(Ef~Ei—uj). 
The multipole coefficients in these sums are the reduced 
matrix elements of the corresponding multipole operators 



CM) = (f\\Cj(q)\\i) 

Lj(q) - (f\\Lj(q)\\i) 

Ej(q) = (fWEjiqW) 

Mj(q) = (/||Mj(g)||i) 



(33) 
(34) 
(35) 
(36) 



The values of J and Jf are related by angular momentum 
conservation \ Ji — Jf\ < J < Ji + Jf ■ In the particular 
case of closed-shell nuclei, such as 12 C, with Ji = 0, we 
have Jf = J. 
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C. Weak charged current 

In order to simplify the comparison with the LFG, in 
this first stage we apply the above formalism to the CC 
= y^—Av j n the static limit. This is not unreasonable 
for the ^-capture reaction since all the momenta involved 
are small. Thus we only maintain the leading order in the 
standard expansion of the matrix clement of the vector 
current 



^(p' > p)=S(p / ) 



2F ^ + iZ 2W° 



2FJ 



«(P) (37) 



in powers of p/M, p'/M, with M the nucleon mass, and 
Q M = (u>, q) the four-momentum transfer (Q 2 = lo 2 — q 2 ). 
Therefore we take 



V° 
V 



2f; 
o 



v 



(38) 
(39) 



In the case of the axial current 

A»(p', p) = u(p') [G A n 5 + GpQ^j 5 ] u(p) (40) 

we expand taking into account that, from PCAC, the 
pseudo-scalar form factor Gp is of order O(M) 



Gt 



2M 



m 2 - Q 



:Ga 



(41) 



and the leading-order term in the expansion of the axial 
current becomes 



.4° 
A 



Ga , v 



ml -Q 2 
G A rr 



G A 



Q 



■(q- <r)q. 



(42) 
(43) 



Therefore the total weak CC in the static limit that we 
use in the present work is 



J 



7° 
•J\- 



V ~ 

J,l 



J° 



p 

Jp 



(44) 
(45) 



and the different terms in these equations are defined be- 
low. First order terms in an expansion in powers of 1 /M 
not included in our calculation can give an appreciable 
contribution, but the present approximation is enough 
for our purposes of testing the equivalence between LFG 
and shell models. 



D. Multipole matrix elements of the current 



The different multipolcs of the vector, axial and 
pseudo-scalar currents (Jy, 3a and Jp), introduced in 
eqs . I|44I45|) . are computed following the approach of ref. 
|34| , where the matrix elements of the electro- weak neu- 
tral current were considered in the context of parity- 
violating electron scattering. In the case of the vector 



current we only consider the zeroth component to lead- 
ing order J v — 2F± . Therefore only the Coulomb multi- 
poles of this current enter in our calculation. The reduced 
matrix elements between single nucleon wave functions, 
with angular momentum quantum numbers (l p ,j P ) and 
(lh, jh), are given by 



(p\\Cj(q)\\h) = 2*f P+ . j[J] ajlj(q). 



(46) 



Here we use the notation _P+ for the parity function (=1 
if n is even and if n is odd), and we have defined the 
function Ij(q) 



dr^jjiq^Rl^Rhir^ir) 



(47) 



which contains the dynamical information on the nuclear 
transition and the muon wave function. Finally the cou- 
pling coefficient aj is defined in terms of a three-j coeffi- 
cient 

flJ = 7^ U -I o J' ( 8) 

In the case of the axial current we only consider the space 
components 3a = G^er (we neglect the time component 
to leading order) so only the longitudinal and transverse 
(electric and magnetic) matrix elements enter: 

(p\\Lj(q)\\h) = iG A P+ +lh+J+1 ^L 

x [(k p + Kh- J)Ij-\{q) 

+ (K P + Kh + J+ 1)7,7+1 (?)] (49) 

(p\\Ef(q)\\h) = -iG A Pt 



ip+ ' h+J+1 y/j(J+l)[J] 

x [(J + 1 + k p + K h )JI J+1 (q) 

+ (J-K p -K h )(J+l)Ij^(q)l50) 

(p\\Mf(q)\\h) = G A P+ +lh+J -0L= 

x (k p - K h )Ij(q), (51) 

where we use the notation k p = (— iyp+ l p+i (j p + i). 
Note that the longitudinal and electric multipoles have 
abnormal parity, i.e., l p + lh + J = odd as expected for 
an axial current. 

In the case of the pseudo-scalar current Jp = 
— m a G J A ga (q • cr)(9 p , the multipoles can be related to 
the longitudinal components of the axial current J Z A — 
G A cr ■ q. Using the expansion (|22l) for the longitudinal 
current we have for the zero-th component 



7° 
■j p 



toq 



Q 



- T z 



bjq 



m 2 - Q 2 



j 



4vr> 'i J [J)L 



A 

JO- 



(52) 
(53) 
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Comparing with the expansion l|21[l we obtain that the 
Coulomb operators of the pseudo-scalar current are pro- 
portional to the longitudinal multipoles of the axial cur- 
rent, 



C p 



JO 



-f A 

)2 JO 



(54) 



and the same relation holds for the matrix elements. 
Since the spatial part of the pseudo-scalar current is pro- 
portional to q, it has no transverse components. Only 
the longitudinal multipoles enter, that are again propor- 
tional to the axial ones 



and cfp = k FP /2M is the Fermi energy of protons, and 
cfn = k FN /2M for neutrons. The LFG width is then 
obtained by inserting the proton and neutron densities, 
Pp(t) and /9/v(r), of the finite size nucleus into eq. (|5fc>|> 
and averaging with the muon density [l| 



Tlfg = J d 3 r\^{r)\ 2 T FG [pp(r),p N (r)}. 



(61) 



An important input for the LFG is the experimental Q- 
value for the reaction 



Q = M{i_ x Y) - M(£X) = LO n 



(62) 



fp _ 



— f A 

)2 JO' 



(55) 



and a similar relation between the corresponding matrix 
elements. 

Finally, note that in the present static approximation, 
where there are no transverse multipoles for the vector 
current, the response function Rt> — 0, because only 
the interference between electric and magnetic multipoles 
of the vector and axial current, respectively, (and vice- 
versa) would enter in eq. J22J). 



E. The local Fermi gas 

In the local Fermi gas model we first compute the decay 
width Tpc[pp, Pn) for a muon at rest inside a Fermi gas 
with constant proton and neutron densities 



which is the minimum value allowed for the energy trans- 
fer lo. In order to account for this value in the Fermi gas, 
we substitute lo by lo — Q, since part of the energy to is 
employed in producing the final nucleus. In this way we 
treat correctly the energy balance, which is important 
for describing the experimental muon-capture width pj. 
When different densities are used for protons and neu- 
trons, especially in the case of 208 Pb, there is a gap 



cfn — £fp 



(63) 



between neutron and proton Fermi energies, that has to 
be considered also in the energy balance by substituting 

lo — > lo + e gap - Q. (64) 



III. RESULTS 



Pp 



"FP 



Pn 



"FN 



/3tt 2 , (56) 



where kpp and k F N are the Fermi momenta of protons 
and neutrons, respectively. With the charged current 
(I44I45|I . the response functions 1|15I|19H arc computed in 
this model using the formalism of [lj, |34| . The final result 
can be written simply as 

" T,<; G2 -e> 2 [AF? v + G 2 A (3 + C P -2Cp)]R (57) 



de' 



7T 



where we have defined the following factor coming from 
the pseudo-scalar current 



C P 



(58) 



and the function R is related to the imaginary part of 
the Linhard function |lll34j 

1 — M 2 
2R Q - — lmU= ——e^FP - e )(e F p ~ e ). (59) 

Here we have defined 



eo = Max < epN — lo, 



1 / Mlo q 
2M 




(60) 



In this section we present results for a set of closed- 
shell nuclei 12 C, 16 0, 40 Ca, and 208 Pb. In the extreme 
shell model the initial and final nuclear wave function 
are described as Slater determinants constructed with 
single-particle wave functions that are solutions of the 
Schrodinger equation with a Woods-Saxon potential 

V{r) = V f(r, Rq, a ) - V LS + V c {r) 



dr 



where 



f(r,R ,a ) = 



1 -|- e (r-Ro)/a 



(65) 
(66) 



and Vc{r) is, for protons, the Coulomb potential of a 
charged sphere of charge Z — 1 and radius Rc, and it is 
equal to zero for neutrons. The parameters of the po- 
tential are commonly fitted to the experimental energies 
of the valence shells or the charge radius. In the present 
case of muon capture we fit the experimental Q-value 
(|62|l for the decay reaction QJ. In the shell model, the 
energy difference between hadronic final and initial states 
is computed as the difference between the corresponding 
shells 



CO 



(67) 
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v p 


Vis 


V N 


Vls 




i2 C WSl 
WS2 
WS3 


-52.38 
-62.38 
-62.38 


-20.30 
-3.20 
-3.20 


-50.85 
-50.85 
-38.30 


-24.11 
-18.40 
-3.15 


1.25 0.57 


ib WSl 
WS2 
WS3 


-52.50 
—52.50 
-50.00 


-0.60 
—7.00 
0.00 


-52.50 
—42.80 
-50.00 


-0.60 
—6.54 
0.00 


1.27 0.53 


4U Ca WSl 
WS2 
WS3 


-50.45 
-57.50 
-57.50 


-4.83 
-11.11 
-11.11 


-48.66 
-55.00 
-53.00 


-5.20 
-2.30 
-5.10 


1.25 0.53 



TABLE I: Parameters of the Woods-Saxon potentials used 
in 12 C, 16 and 40 Ca for protons (P) and neutrons (N). The 
units are MeV for V;, and fm for ao and ro. The reduced 
radius parameter ro is defined by Rq = roA 1 ^' 1 . The Coulomb 
radius is chosen Rc = Ro- 







V 


V LS 


ro 


ao 


tls 


aLS 


WSl 


P 


-60.4 


-7.45 


1.26 


0.79 


1.21 


0.59 




N 


-46.9 


-5.64 


1.21 


0.66 


1.17 


0.64 


WS2 


P 


-60.4 


-6.75 


1.26 


0.79 


1.22 


0.59 




N 


-43.5 


-6.08 


1.26 


0.66 


1.17 


0.64 



TABLE II: Parameters of the Woods-Saxon potentials of 
208 Pb for protons (P) and neutrons (N). Note that we use dif- 
ferent radius parameters for the central and spin-orbit parts of 
the potential. The units are MeV for Vi, and fm for a^ and Vi. 
The reduced radius parameters Ti are defined by Ri = TiA 1 ^ 3 . 
The Coulomb radius is chosen Rc = Ro- 

where e p and eh are eigenvalues of the Schrodinger equa- 
tion for particles (neutrons) and holes (protons), respec- 
tively. Therefore the Q value I)(i2ll is obtained in this 
model as the energy difference between the first unoc- 
cupied neutron shell and the last occupied proton shell, 
corresponding to the transition of a valence proton to a 
neutron above the Fermi level. This makes only one con- 
dition for fixing the several parameters of the potential 
(|65[) . Wherever possible we set the remaining parameters 
of the potential to values similar to the ones used in other 
studies like those of refs. pH 133. l35j | . In our calculation 
we use different sets of parameters, denoted WSl, WS2 
and WS3, shown in table Qj for 12 C, 16 0, 40 Ca, and in 
table El for 208 Pb. 

The only states relevant for /i-capture are the occupied 
proton holes and the neutron particles above the valence 
shell. In the discrete sector several transitions are pos- 
sible with fixed excitation energies. The single particle 
energies of the last occupied proton shell and first un- 
occupied neutron shell obtained with the potentials of 
tables IIIIII are shown in table IIIII The Q value corre- 
sponds to the transition P — ► N in table 11111 with an 
energy difference 

Q = e(N) - e(P) (68) 

which is also shown in table IIIII together with the ex- 
perimental value Qoxp in the last column. The number 
of discrete neutron states is finite. Above the last dis- 
crete neutron state, the next allowed transitions are to 



Nucleus WSl WS2 WS3 esp 



1Z C 


■P1P3/2 
iVlpi/2 

Q-value 


-15.96 
-2.08 
13.88 


-18.38 
-4.50 
13.88 


-18.13 
-4.25 
13.88 


13.880 


IbQ 


Plpi/2 

Nld 5/2 
Q-value 


-15.31 
-4.39 
10.92 


-12.77 
-1.84 
10.93 


-13.76 
-2.83 
10.93 


10.931 


La 


Nlf V 2 
Q-value 


— O.OO 

-6.51 
1.83 


—a. to 

-6.95 
1.83 


— O. (O 

-6.95 
1.83 


1.822 


2U8 Pb 


P3s 1/2 
N2g 9/2 
Q-value 


-8.19 
-2.68 
5.51 


-8.19 
-2.68 
5.51 




5.512 



TABLE III: Single particle energies in MeV used in the fit of 
the Q-value for ^-capture (the experimental values are shown 
in the last column). 







discrete 


total 


LFG 


% 


12 C 


WSl 


0.3115 


0.4406 


0.4548 


3.2 




WS2 


0.3179 


0.4289 


0.4360 


1.7 




WS3 


0.2746 


0.5510 


0.4732 


-14.1 




WSl 


1.113 


1.282 


1.360 


6.1 




WS2 


0.590 


1.118 


1.392 


24.3 




WS3 


1.154 


1.332 


1.387 


4.1 


4U Ca 


WSl 


29.10 


37.12 


36.73 


-1.1 




WS2 


27.79 


33.79 


34.90 


3.3 




WS3 


26.28 


32.73 


35.03 


7.0 


2U8 Pb WSl 


215.6 


390.3 


399.4 


2.3 




WS2 


266.8 


467.4 


439.5 


-5.9 



TABLE IV: Integrated width in units of 10 5 s _1 for the dif- 
ferent nuclei and Woods-Saxon potentials, compared with the 
LFG results using the corresponding charge densities. The 
discrete contribution of the shell model is shown in the first 
column. 

the continuum. The continuum neutron states are ob- 
tained by solving the Schrodinger equation for positive 
energies. More details on the continuum solutions can 
be found in refs. 0111. 

In order to compare with the LFG, it is important to 
use as input the proton and neutron densities obtained in 
the corresponding shell model, by summing over occupied 
states as follows: 

P p ^= E ^-\Rn h (r)\ 2 , (69) 

protons 

where R n ij (r) are the radial wave functions, and a similar 
expression for neutrons. 

The shell model calculation has been checked by com- 
parison with the factorized plane wave impulse approxi- 
mation (PWIA) |36j |. In this approximation there is no 
final-state interaction and hence the final neutron states 
are plane waves. The transition matrix elements appear- 
ing in the hadronic tensor l|13f) are computed trivially in 
terms of the product of a single nucleon current matrix 
element times the Fourier transform of a nuclear over- 
lap function of the missing momentum, in a way which 
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wsi 

LFG 



FIG. 1: Proton and neutron densities of C for the several 
WS potentials used in this work. 



is similar to the analysis of exclusive (e, e'p) reactions 
|37| (but this time the nuclear overlap function includes 
the bound muon wave function). As a consequence, the 
exclusive hadronic tensor factorizes as the product of a 
single nucleon hadronic tensor times a partial momentum 
distribution, and the calculation is straightforward in the 
shell model. For the present case an additional integra- 
tion and a sum over initial states is needed since we are 
interested in the inclusive case, similar to the factorized 
PWIA in (e,e') introduced in [sjj. The PWIA can be 
also approached with our multipole expansion code by 
setting to zero the WS potential in the final states. This 
allows us to check the multipole expansion calculation 
and, at the same time, to fix the number of multipoles in 
the sum over J, eqs. <|2 II 12311 . The differences with the 
factorized calculations are negligible when we include up 
to five multipoles. 

In table Hvl we show results for the integrated inclusive 
widths for the four nuclei and for the different models 
used in this work. For each one of the WS parameteri- 
zations we show in the second column the contribution 
to the width from the discrete final neutron states, while 
in the third column we show the total width (discrete + 
continuum). The LFG results are shown in the fourth 
column, and for comparison we show the percentual rel- 
ative difference between LFG and WS in the last column. 
Next we discuss the results obtained for each one of the 
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FIG. 2: Differential SM width of 12 C to the continuum (left 
panels) compared to the LFG, and partial widths to the dis- 
crete states (right panels), as a function of the neutrino en- 
ergy, for the different WS potentials considered in this work. 



nuclei studied in this work. 



In tabic II VI we can see that, in the case of WSI and 
WS2, the LFG and WS results for 12 C are quite simi- 
lar, differing only in ~ 2-3%. In the case of WS3 the 
differences are larger, around 14%. 

All the WS potential parameters have been fitted to 
the experimental Q- value (|68|) . which, in the SM, is the 
difference between the neutron P1/2 and proton p 3 / 2 ener- 
gies. Among all the potentials, WSI is the more realistic 
since it also fits the proton and neutron separation en- 
ergies of 12 C and give reasonable masses for the ground 
states of the 13 N and 13 C nuclei. When we use similar pa- 
rameters for protons and neutrons, like in WSI, we need 
a large spin-orbit splitting in order to fit the experimental 
Q- value. In the case of the potential WS2 we use different 
parameters for protons and neutrons: The proton well is 
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FIG. 3: The same as Fig. 1 for 10 O. 



similar to the one of ref. [3j|, that is more attractive 
than WS1, with small spin orbit strength. The neutron 
parameters are similar to WS1. Finally in WS3 we have 
used a small neutron spin-orbit splitting, as for protons, 
but we had to make the neutron well much less attractive 
than for protons. Apart from changing the single-particle 
energies, the effect of modifying the WS potential can be 
appreciated in the proton and neutron densities shown 
in Fig. 1. For more attractive potentials the nucleus be- 
comes more dense in the interior. For this reason, the 
WS3 neutron density turns out to be the smallest one, 
while the proton density is around 3/2 the neutron one. 
Hence the LFG results are worse for very different neu- 
tron and proton densities. In this situation the proton 
and neutron Fermi momenta (|56[) are clearly different, 
leading to a gap between proton and neutron energies, 
(|63|1 . which in this case is negative, since the density is 
smaller for neutrons, and €fn < epp- Therefore a proton 
near the Fermi surface can decay to a neutron above the 
neutron Fermi surface with an energy decrement. This 
is an unrealistic situation, since precisely in this case the 
neutrons are less bound than protons in the SM, and 
therefore lie at higher energies. Another argument to 
disregard this case is the well known property of closed 
(sub)shcll nuclei such as 12 C, for which the neutron and 
proton densities should be similar. Note that in all cases 
the gap between the N and P Fermi species has been 
taken into account in the energy balance by the replace- 
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ment (lf)4Tl . 

In Fig. 2 we compare the SM results for the differential 
width to the continuum with the LFG distribution for 
the different WS parameters (left panels). The shapes 
of both distributions are completely different. The dif- 
ferences are more apparent for the WS1 potential, where 
there is a very high and sharp neutron resonance in the 
SM spectrum. The partial widths to the discrete states 
are shown in the right panels of Fig. 2. Considering these 
differences in shape between the LFG and the SM, it is 
a very notable result that the integrated widths (adding 
the discrete states) take similar values in both models as 
was shown in table IIVI This outcome agrees with the 
findings of ref. [38j, where the same problem was ad- 
dressed in the context of inelastic electron scattering on 
nuclei. 

The biggest contribution to the width comes in all the 
cases from the transition to the ground state, and its 
magnitude does not depend very much on the potential, 
since in the transition p 3 / 2 — > P1/2 the wave functions 
in the initial and final states are similar across the dif- 
ferent potentials. Note also that there are transitions 
to final discrete states that lie in the continuum (par- 




FIG. 5: The same as Fig. 1 for 40 Ca. 



ticularly, transitions from the Is shell). These states 
will contribute to the giant resonances after an appropri- 
ate treatment of the residual interaction (such as in the 
RPA). Under the light of the present preliminary study 
and the results of ref. pj one expects that the inclusion 
of the RPA does not change too much our conclusions 
and the total integrated width be similar for correlated 
LFG and SM. 



O 



In the case of 16 O the integrated widths computed in 
the LFG are also very close, ~ 4-6%, to the SM results 
with the potentials WS1 and WS3 (see table E2J- The 
worse results are obtained for the WS2 parameterization; 
the corresponding width is 24% of the SM one. This can 
also been understood in terms of what was said for the 
case of 12 C above, by looking at the 16 densities shown 
in Fig. 3. The case of WS2 is the only one where the 
protons are more bound than the neutrons, hence the TV- 
density is smaller than the P-density inside the nucleus, 
which is again an unrealistic situation because one expect 
the opposite in a closed shell nucleus such as 16 O. 

We should add that the 16 O nucleus is delicate in the 
sense that the experimental Q-value of 10.93 MeV is too 
large to be fitted by the WS parameters found in the liter- 
ature 0, 01 • I n fact m the SM the Q- value is the differ- 
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FIG. 6: The same as Fig. 2 for 40 Ca. 



ence between the and -Plpi/2 energies (see table 

llllf) . The effect of the spin-orbit potential is to increase 
e (Pi/2) and to decrease e(<2 s / 2 ), that is, goes to reduce 
the Q-value. (The opposite happens for 12 C, where the 
Q- value is the difference between the Npi/2 and Ppz/i 
energies. Hence the spin-orbit goes to increase the Q- 
value.) Therefore to make that value as large as 11 MeV 
one needs a small spin-orbit potential, as in WSl, or to 
rise the neutron well with respect to the proton well, as 
in WS2, at the cost of making the neutrons less bound 
than protons. The first option is preferred because it al- 
lows for similar proton and neutron densities. Precisely 
the third parameterization WS3 has been chosen with 
Vls — to maximize the difference between these two 
shells. 

The differential and partial widths of 16 O for the three 
WS potentials are shown in Fig. 4. In the three cases the 
main contribution comes from the discrete spectrum (see 
also table IIV(I . Since Vls — for WS3, the dominant 
contribution comes from transitions from the lp to the 
Id shell. 
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The LFG results improve when the mass of the nucleus 
increases as in the present case of the nucleus 40 Ca. In 
fact, from table II V I we see that for this nucleus the LFG 
integrated width is within 1% of the SM result for WS1, 
and 3% and 7% in the other two cases. This improvement 
was expected because the Fermi gas description of the 
nucleus should work better for heavier nuclei. In the case 
of WS2 and WS3 the proton parameters have been fixed 
to the typical values used in the literature, and we have 
fitted the neutron ones. Since here the experimental Q- 
value is small, Q — 1.8 MeV, one does not need to change 
too much the typical neutron parameters. In the case of 
WS1 we have tried to maintain the P and N parameters 
similar. The proton and neutron densities are close in all 
cases, as shown in Fig. 5, and the proton levels always 
lie above the neutron ones. 

The neutrino spectrum shown in Fig. 6 presents a 
more complex structure than the lighter nuclei discussed 
above. More potential resonances arise and the discrete 
spectrum presents more lines distributed along the al- 
lowed energy region. 

D. 208 Pb 

Finally we discuss the results for the closed-shell heavy 
nucleus 208 Pb. In table II Vl we present integrated widths 
only for two sets of potential parameters, WS1 and WS2. 
This is the only case where we use different radius param- 
eters for protons and neutrons and also for the central 
and spin-orbit parts of the potential, see table [HI In 
both cases the LFG results are close, within 3 and 6%, 
to the SM ones. The Q value, 5.5 MeV, is close to the 
N2g Q /2 and P3s 1/2 energy difference of typical param- 



FIG. 8: The same as Fig. 2 for 208 Pb. In this case we show 
the discrete contributions in units of 10 4 s _1 in the same panel 
as the continuum one. 

eterizations [3 HH • Only small variations of these pa- 
rameterizations found in the literature are allowed if one 
wants to maintain the ordering of the energies around 
the Fermi level. Also only small variations are needed 
to fit the experimental Q-value. In the present case the 
treatment of asymmetric nuclear mater is essential, be- 
cause the proton and neutron densities, shown in Fig. 7, 
are clearly different. Therefore the correct treatment of 
the gap in the energy balance, eqs. I|63I64|) . is needed 
to obtain the results of table llVl Moreover in this case 
the neutrino spectrum shown in Fig. 8 shows also an im- 
proved resemblance between the LFG and SM (although 
numerous potential resonances appear), even taking into 
account the distribution of the discrete spectrum. 

IV. CONCLUSIONS 

In this paper we have estimated the magnitude of the 
finite nucleus effects on inclusive muon capture, aiming at 
quantifying the uncertainty of the LFG results of ref. [l| . 
It is not possible to disentangle these effects by compar- 
ison with the highly sophisticated RPA or shell models 
existing in the literature due to the different theoretical 
ingredients embodied in them. 

To know how much the LFG is modified by finite size 
effects, one would need a finite nucleus model with ex- 
actly the same input as the LFG, in order to make the 
comparison meaningful. Obviously this would be a draw- 
back because precisely one wants to use the LFG due to 
its simplicity, in order to include very complex dynamical 
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effects hard to incorporate in finite nucleus treatments. 
Therefore before using a very sophisticated model, it is 
convenient to see what happens in the uncorrelated case. 

In this paper we have focused on a simple shell model 
without nuclear correlations, but that contains the rele- 
vant information about the finite nuclear structure, and 
we have compared it with the uncorrelated LFG using 
the same input. In particular the SM proton and neu- 
tron densities have been used in the LFG calculation. We 
have applied both models to a set of closed-shell nuclei: 
i2 Cj 16 0) 40 C&j and 208p k In the SM we fit the experi- 
mental Q-value of the decay, while the same value is used 
to correct the energy transfer in the LFG, taking into ac- 
count also the gap between neutron and proton Fermi 
energies. As expected, the neutrino spectrum is very 
different in both models, in particular the LFG cannot 
account for the resonances and discrete states. However, 
in the case of the lighter nuclei, 12 C and 16 0, the SM and 
LFG results for the integrated width are close — within 3— 
6% — for WS parameters with similar neutron and proton 



densities, but the results are somewhat different, within 
14-24%, for the disregarded cases in which the protons 
lie below the neutrons. For the medium and heavy nu- 
clei, 40 Ca and 208 Pb, the integrated widths are always 
very close, within 1-7%. The final neutrino spectra of 
the LFG becomes more similar to the SM, including the 
discrete part, for heavier nuclei. Under the assumption 
that RPA correlations and finite size effects are some- 
what decoupled for integrated inclusive observables, the 
present results can explain why the LFG results of ref. 
PJ describe so well the experimental data. 
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